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Abstract 

The effect of CP violating phases in the MSSM on the relic density of the lightest super- 
symmetric particle (LSP) is considered. In particular, the upper limits on the LSP mass are 
relaxed when phases in the MSSM are allowed to take non-zero values and when the LSP is 
predominantly a gaugino (bino). Previous limits of < 250 GeV for Qh 2 < 0.25 can be 
relaxed to mg < 650 GeV. The additional constraints imposed by the neutron and electron 
electric dipole moments induced by these phases are also considered. Though there is some 
restriction on the phases, the bino mass may still be as large as ~ 350 GeV and certain 
phases can be arbitrarily large. 
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It is well known that by considering the cosmological relic density of stable particles, one 
can establish mass limits on these particles. In the minimal super symmetric standard model 
(MSSM), because of the many unknown parameters it is a significantly complicated task to 
set limits on the mass of the lightest particle, the LSP. The annihilation cross section will 
depend on several parameters which determine the identity of the LSP as well as parameters 
such as scalar masses which mediate LSP annihilations. In addition, there are (at least) two 
CP violating phases in the MSSM which play a role in determining the LSP relic abundance 
[[[]]. As will be shown, when the phases are allowed to take their maximal values, mass 
limits on the LSP are greatly relaxed. However, constraints from the neutron electric dipole 
moment will restrict the degree to which these mass limits can be relaxed. 

The mass and relic density of the LSP depends on a relatively large number of parameters 
in the MSSM. With some assumptions on these parameters, their number can be whittled 
down to 6: 1) if the gaugino masses are assumed to obey a GUT relation, they can be 
characterized by a single mass parameter, M 2 ; 2) the Higgs mixing mass, /i; 3) the ratio of 
Higgs vevs, tan/3; 4) the soft (diagonal) sfermion mass parameters, mj (assumed here to be 
degenerate at the weak scale); 5) the soft trilinear terms, Af (also assumed equal at the weak 
scale); and 6) the soft bilinear, B. In general, all but the sfermion masses, mj, are complex. 
However not all of the phases associated with these are parameters are physical 0. It is 
common to rotate away the phase of the gaugino masses, and to make -B/i real, which ensures 
that the vacuum expectation values of the Higgs fields are real. For simplicity, generation 
mixing in the sfermion sector will be ignored though left-right mixing will be included. If 
furthermore we assume that all of the A/'s are equal, and simply label them by A, we are 
left with two independent phases, the phase of A, 9a, and the phase of /x, 9^. The phase of 
B is fixed by the condition that _B/i is real. 

There has been a considerable amount of work concerning phases in the MSSM. For 
the most part these phases are ignored because they tend to induce large electric dipole 
moments for the neutron |2|. If one considers the direct contribution to the edm from 
gluino exchange, one finds a dipole moment 

Q. s Tfl u 7Tlq \Tfl t \ . 1 n— 23 • f -t \ 

a n ~ e -. sm7 t ~ few x 1(J sin 74 ecm (1) 

where the gluino mass, m- g and mj as well as m t = A* + /icot (3 are all taken to be of order 
100 GeV. The angle 7^ = argm t is clearly a function of both 9^ and 9a- To suppress the 
electric dipole moments, either large scalar masses (approaching 1 TeV) or small angles (of 
order 10~ 3 , when all SUSY masses are of order 100 GeV) are required. For the most part, 
the community has opted for the latter, though the possibilities for large phases was recently 
considered in M. To reconcile large phases with small electric dipole moments, some of the 
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sparticle masses are required to be heavy. In ||, it was first assumed that the SUSY phases 
were large (7r/4), and by computing the various contributions to the edm of the neutron from 
neutralino, gluino, and chargino exchange (it was the latter that they found to be dominant) 
either large sfermion or neutralino masses (or both) were required. For example, Kizukuri 
and Oshimo Q found that for M 2 — fi — 0.5 TeV, and tan/3 = 2 the chargino contribution 
to the edm of the neutron exceeds the experimental upper bound || of \d n \ < 1.1 x 10 25 ecm 
unless, rrij > 2 TeV. When M 2 = \x = 0.2 Tev the limit on the sfermion masses becomes 
in? > 3 TeV for tan ft = 2 and m s > 6 TeV for tan ft = 10. 



Figure l:Typical contribution to the edm. 

However, unless .R-parity is broken and the LSP is not stable, one requires that sfermions 
be heavier than neutralinos, and if they are much heavier, this would result in an excessive 
relic density of neutralinos. Thus one must determine the relationship between potentially 
large phases in the MSSM and the relic density while remaining consistent with experimental 
bounds on the electric dipole moments. Here I will consider only neutralinos as the LSP. 
The mass limits will then depend sensitively on the masses of the scalars. For example when 
parameters are chosen so that the LSP is a gaugino (when the supersymmetry breaking 
gaugino masses Mi$ are taken to be smaller than the supersymmetric Higgs mixing mass, 
/i), for a given LSP mass, the requirement that Qh 2 < 1/4 places an upper bound on the 
sfermion mass (Q is the fraction of critical density and h is the hubble parameter scaled to 
100 km/s/Mpc). The reason being that flh 2 oc 1/ (<jv)a, and {<tv)a oc mj 4 is the annihilation 
cross section. There is of course in addition, a dependence on the gaugino mass such that 
for m gaugino < rrif, a larger gaugino mass implies a larger annihilation cross section. 

For small gaugino masses, the LSP is a photino, and one obtains a lower limit on the 
photino mass |]]. At larger gaugino masses the LSP is a bino whenever M 2 < \x and 
M 2 > 200 GeV; for smaller M 2 , the bino is still the LSP for large enough values of //. The 
bino portion of the M 2 — \x LSP parameter plane is attractive, as it offers the largest possibility 
for a significant relic density |7], || 0- The complementary portion of the parameter plane, 
with fi < M 2 , only gives a sizable density in a limited region, due to the large annihilation 
cross sections to W + W~ and ZZ and due to co- annihilations [1(J with the next lightest 



neutralino (also a Higgsino in this case), which is nearly degenerate with the LSP |l.l|. For 



the remainder of the discussion here, I will focus on binos as the LSP. In addition to resulting 
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in a sizable relic density, the analysis is simplified by the fact that in the nearly pure bino 
region, the composition and mass of the LSP is not very sensitive to the new phases. However, 
as is shown, the relic density, which is determined primarily by annihilations mediated by 
sfermion exchange, is quite sensitive to the phases, j t and 7&. 

Typically, bino annihilation is dominated by the BB — > ff channel. When the thermally 



averaged cross-section is expanded in powers of T/m-g (see e.g. |jl2|), the cross section can 
be written as 

{av) a = a + bT / nig (2) 

where 



a oc 



— (Yl + Y*Y^( (3) 
128vr \ L R J Aj K J 

b-^iYt + Y^)^ (4) 
128vr y L R J A} K ' 

and A 2 = mj + m 2 ~ — m 2 . In the above equations, Yl(b) is the hypercharge of fi(R). The 
expression for b is significantly more complicated but what is important to note is that b 
contains a piece which is proportional to m 2 ~ whereas a does not. This is a manifestation of 
the p-wave suppression of the annihilation cross-section for majorana fermions. 

When common sfermion masses are assumed and when sfermion (left-right) mixing is 
ignored, there is an upper limit < 250-300 GeV [0, in order to ensure that Q B h 2 < 1/4. 
(This upper limit is somewhat dependent on the value of of the top quark mass. In |7|, [13[] , 
we found that the upper limit was ~ 250 GeV for m t ~ 100 GeV. When m t = 174 GeV, 
the upper limit is 260 GeV. For m t ~ 200 GeV, this limit is increased to < 300 GeV. 
Furthermore, one should be aware that there is an upward correction of about 15% when 



three-body final states are included |14[ which raises the bino mass limit to about 350 GeV. 
This latter correction would apply to the limits discussed below though it has not been 
included.) As the bino mass is increased, the sfermion masses, which must be larger than 
mg, are also increased, resulting in a smaller annihilation cross section and thus a higher 
relic density. At ~ 250 GeV, even when the sfermion masses are equal to the bino mass, 
flh 2 ~ 1/4. Note that this provides an upper bound on the sfermion masses as well, since 
the mass of the lightest sfermion is equal to the mass of the bino, when the bino mass takes 
its maximum value. One should also note that in the range 20 GeV < < 300 GeV, it is 
always possible to adjust nij in order to obtain Q B h 2 = 1/4. 

In [[13], it was shown that the upper limit to the bino mass is sensitive to the level of 
sfermion mixing. The general form of the sfermion mass 2 matrix is |TH 



Ml + m 2 + cos2/3(T 3/ - Q f sin 2 6 W )M 2 m f |r%|e^ 

m f \m f \e-^f M 2 R + mj + cos 2(3Q } sin 2 9 W Mf . 



(5) 
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where Ml(r) are the soft supersymmetry breaking sfermion mass which are assumed to be 
generation independent and generation diagonal and hence real. I also assume Ml = Mr. As 
discussed above, there is a non-trivial phase associated with the off-diagonal entries, which 
is denoted by m f (\m f \e' l ' yf ), of the sfermion mass 2 matrix, and 

\Mf\e*" = Rffi + A* = Rf\/J,\e i9 » + \A\e~ l9A (6) 

where Rf = cot/3 (tan/?) for weak isospin +1/2 (-1/2) fermions. 



Figure 2: Upper and lower limits to the bino mass as a function of Af for \i = —3000 GeV. 



In general, one finds for the s-wave annihilation cross-section |l ] 



a f 



1287T 



Ai(rrifW 2 — 2rrigWZ + rrifZ 2 ) + A2(m/x 2 + 2m^xy + rrify 



A X A 2 



(7) 



where A» = m 2 ? + m~ — m 2 and the factors x, y, w, z are determined from the bino-fermion- 



B 



sfermion interaction Lagrangian, 



C ffB = T2 9 ' (YRfPLBf R + Y L fP R Bf L ) +h.c. 



= j^g'f [fixP L + f2wP L + hyP R - f 2 zP R ) B + h.c. (8) 

When phases are ignored, x = Y R sm9f, y = Y L cos9f, w = Y R cos9f, z = Y L sm9f. 
The angle 9j relates the sfermion mass eigenstates to the interaction eigenstates by f\ = 
Jl cos 9 f + f R sin 9 f and f 2 = - f L sin 9 f + f L cos9 f . 

When the off diagonal elements of the sfermion mass matrix vanish = 0), Ai = A2 

and xy = wz and the pieces proportional to in (0) cancel yielding the suppressed form of 
the s-wave cross-section given in @. When sfermion mixing is included [0], the limits, which 
now depend on the magnitude of the off-diagonal elements mffrif, are modified. In Figure 
2, the upper and lower limits to the bino mass is shown as a function of Af, which for fixed 
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fi is equivalent to showing the limits as a function of the off diagonal element \m\. Because 
it is the combination mfm s that appears in the off diagonal element, sfermion mixing really 
only affects the stops. Thus when mixing is important one of the stop masses is much less 
than the diagonal elements M L ^ R y In order to ensure that the bino is lighter than the light 
stop, the diagonal elements must be pushed up relative to the no mixing case. This results 
in a reduction in the annihilation cross-section away from the no stop mixing region near 
A = 1500 GeV. The upper limit is asymmetric about this point as the couplings of the bino 
to t\ and t 2 differ, and the lightest stop is either t\ or t 2 , depending on whether or not A is 
greater than or less than 1500 GeV. The lower limit on is easily understood: when there 
is a lot of stop mixing and the diagonal elements are large, the annihilation into channels 
other than tops, is suppressed. If mg < m t , the annihilation cross-section is too small and 
Qh 2 is too large. Thus we must require the top channel to be available and mg > m t . Near 
the no stop mixing region, we have assumed that > 74 GeV |16 |. In addition, we are able 



to obtain Qh 2 = 1/4 for a wide range of bino masses (as was the case when sfermion mixing 
was ignored by varying the the sfermion mass) as well as for a wide range of sfermion masses, 
by varying the magnitude of the off diagonal elements in the sfermion mass matrix. For fixed 
\i this corresponds to varying A. The upper limit to the bino mass is relaxed considerably 
when the phases are allowed to take non-zero values. 

As we have just seen that though sfermion mixing has a large effect on the bino mass 
limits, the mass limit was not greatly increased. Because of the mass splitting among the 
sfermions, (A 2 Ai) there is a reduction in the cross-section which compensates somewhat 
for the presence of the term proportional to in a. With non-zero phases, the enhancement 
in the cross-section is much greater. Even with A 2 ~ Ax, the terms of interest in (0) combine 
to 2m^(xy — wz) ^ which is in fact proportional to shi7y. The effect of the phases on 
the bino mass limits is shown in Figures 3, 4 as a function of the magnitude of the off- 
diagonal term in the top-squark mass 2 matrix, m t , given the conditions: 1) Q^Ii 2 < 1/4, 
2) the lightest sfermion is heavier than the bino, and 3) the lightest sfermion is heavier 
than 74 GeV. In both figures we have taken tan (3 = 2 and m top = 174 GeV. In Figure 3, 

= 3000 GeV and in Figure 4, \fi\ = 1000 GeV. The various curves are labeled by the 
value of 7;, assumed, and in addition, has been maximized for all allowed values of 9^. 
The lower limit on assumes 75 = n/2. As one can see, when 7^ is allowed to take its 
maximal value of tt/2, the upper limits are greatly relaxed to < 650 GeV. With and 
7b fixed, for a given value of m t and 9^ all of the remaining quantities such as \A\, 9a, j t , and 
m b are determined. 

For large m t , the diagonal mass terms M| must be taken large to ensure that the mass of 
the lightest stop is > mg. This drives up the masses of the other sfermions and suppresses 
their contribution to the annihilation. As m t is decreased, M| must drop and the other 
sfermions begin to contribute and the upper bound on is increased. In particular, 
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annihilation to /i's and r's becomes important, since 



c,t 



YlYl 



s,6 



\\ : 4 : 1 



(9) 



Decreasing 7 & reduces the effect of /i's and r's, and this can be seen as a decrease in the upper 
bounds in Figures 3 and 4. For fn t sufficiently small, the stops become unmixed, diminishing 
somewhat their contribution and slightly decreasing the upper bound on mg. 



Figure 3: Upper limits on the bino mass as a function of the off- diagonal element m t in the top squark 
mass 2 matrix, for various values of 75, the argument of the off-diagonal element of the T3 = — 1/2 
sfermion mass 2 matrix. Also shown is the lower bound (lowest curve) on the bino mass assuming 
7b = 7r/2. The value of |/i| was chosen to be 3000 GeV. 



Figure 4: As in Fig. 3, with |/i| = 1000 GeV. 

I turn now to the calculation of the electric dipole moments of the neutron and the 
electron. As discussed above, the edm's of the electron and quarks receive contributions from 
one-loop diagrams involving the exchange of sfermions and either neutralinos, charginos, or 
(for the quarks) gluinos. In the case of the neutron edm, there are additional operators 
besides the quark electric dipole operator, 7 = ^qa flu qF' lu |3j which contribute. They 
are the gluonic operator Oq = —hf abc G a GbG c |HJ and the quark color dipole operator, 
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O q = jqa^ u qT a G^ u |Tj|. The gluonic operator is the smallest [18] when all mass scales 
are taken to be equal. These three operators are conveniently compared to one another in 
PU[ and relative to the gluino exchange contribution to the 7 operator, it was found that 
7 : O q : O g = 21 : 4.5 : 1. 

Because of the reduced importance of the additional operators contributing to the neu- 
tron edm, it is sufficient to only include the three contributions to the quark electric dipole 
moment. The necessary CP violation in these contributions comes from either in the 
sfermion mass matrices or # M in the neutralino and chargino mass matrices. Full expressions 
for the chargino, neutralino and gluino exchange contributions are found in ||. The depen- 
dencies of the various contributions on the CP violating phases can be neatly summarized: 
for the chargino contribution 

d^-sin^, (10) 
with essentially no dependence on 7/; whilst for the gluino contribution 

d'j ~ m f shi7j , (11) 

independent of and the neutralino contribution has pieces that depend on both sin # M and 
7% sin 7/. All three contributions can be important (including the neutralino contribution, 
in the case of the electron edm), and depending on sin^ and sin7j, they can come in with 
either the same or opposite signs. In particular, signfrfj /d^] = signfsin^/ sin 7/]. For the 
mass ranges we consider, the dipole moments fall as the sfermion masses are increased, and 
sfermion masses in the TeV range can bring these contributions to the neutron and electron 
electric dipole moments below the experimental bounds of |<i n | < 1.1 x 10 _25 e cm and 
\d e \ < 1.9 x 10~ 26 e cm |2~I|, even for large values of the CP violating phases M. However, 



these large sfermion masses are inconsistent with the cosmological bounds mentioned above, 
where sfermion masses must be relatively close to the bino mass in order to keep the relic 
density in check. 

To determine the allowed parameter ranges [l|], we first fix the value of 7^ and take 
\fx\ = 3000 GeV. Then for several values of m t between and 1500 GeV, we determine the 
upper bound on mg, as a function of 6^. As we vary # M across its full range, m b and 7t 
change, and this affects the annihilation rate and consequently the bound on mg. Taking 
at its maximum value allows us to take M\ as large as possible; although the electric 
dipole moments depend on as well, the dependence on M| is sufficiently strong that 
the edm's take their minimum values for the maximum values of and M£. The quark 
and electron edm's can then be computed as a function of 9^ and m t , and one can use the 
nonrelativistic quark model to relate the neutron edm to the up and down-quark edm's via 

d n = ^(4d d -d u ). (12) 

If there is no region of the 6^-m t parameter space which satisfies both the neutron and 
electron edm bounds, we can decrease 7& and repeat the procedure. In practice, the bound 
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on the neutron edm the more difficult of the two to satisfy, and every region of the parameter 
space we show which produces an acceptable neutron edm also produces a sufficiently small 
electron edm. 

For the large value of = 3000 GeV, the largest contribution to the neutron edm comes 
either from gluino exchange (for the more negative values of 9^) or chargino exchange (for 
the more positive values of 6^), and the value for \d n \ is too large unless % takes a relatively 
small value. In particular, non-negligible experimentally acceptable regions of the parameter 
space are found only for 75 < n/25. In Figure 5, a contour plot of the neutron edm as a 
function of 9^ and m t for 7^ = 7r/40 is shown |l|]. The shaded regions demarcate the range of 
# M for this choice of 7&. Much of this range produces a sufficiently small \d n \. As we increase 
m t , the d and u masses become large and \d n \ falls. As we move to values of m t greater than 
~ 1500 GeV, we begin to require a significant tuning of M\ to produce Q^h 2 < 1/4. 



Figure 5: Contours of the neutron electric dipole moment, d n , in the 9^ - WE t plane in units of 10~ e cm. 
The value of was chosen to be 3000 GeV and ■yj, = 7r/40. The shaded region corresponds to values of 
9^ and fn t which are not allowed algebraically for this value of \X and 7fe. 

This procedure can be repeated for \\x\ = 1000 GeV and the result is shown in Figure 
6 0. For lower \fi\, the chargino exchange contribution is enhanced relative to the gluino 
exchange contribution. In this case, one can allow larger values of 7^ up to about tt/6. For 
the same range in m t , 9^ can takes values from -0.3 to 0.2. For either value of /z, the angles 
7t and 9a are unconstrained. 

In summary, we have found that CP violating phases in the MSSM can significantly 
affect the cosmological upper bound on the mass of an LSP bino. In particular, taking the 
maximal value ir/2 for the phase 76 of the off-diagonal component of the T3 = —1/2 sfermion 
mass matrices pushes the upper bound on mg up from ~ 250 GeV to ~ 650 GeV. When we 
additionally consider constraints on neutron and electron electric dipole moments, we find 
the upper bound on m~ is reduced to ~ 350 GeV. Various combinations of the CP violating 
phases are constrained as well: \9^\ < 0.3 and |7&| < n/6 for \fi\ > 1000 GeV, while 74 and 
9a are essentially unconstrained. We note that although the bounds on # M and 7^ are small, 
they are much larger than the values of order 10~ 3 typically considered. 
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Figure 6: As in Figure 5, with chosen to be 1000 GeV and 7& = 7T/8. 
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